Mirror matter can alleviate the cosmological lithium problem 
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The abundance of lithium-7 confronts cosmology with a long lasting problem between the pre- 
dictions of standard big bang nucleosynthesis and the baryonic density determined from the cosmic 
microwave background observations. This article investigates the influence of the existence of a 
mirror world, focusing on models in which neutrons can oscillate into mirror neutrons. Such a 
mechanism allows for an effective late time neutron injection, which induces an increase of the de- 
struction of beryllium-7, due to an increase of the neutron capture, and then a decrease of the final 
lithium-7 abundance. Big bang nucleosynthesis sets constraints on the oscillation time between the 
two types of neutron and the possibility for such a mechanism to solve, or alleviate, the lithium 
problem is emphasized. 

PACS numbers: 98.80.-k, 26.40.-|-r, 14.20.Dh, 12.60.-i 



I. INTRODUCTION 

The abundances of the light elements produced during 
the primordial nucleosynthesis (BBN) in the early hot 
phase of the Universe is one of the historical pillar of the 
big-bang model 1] . The prediction of these abundances 
depends on: 

1. the theory describing gravitation, which determines 
the cosmic expansion rate during BBN via the 
Friedmann equation, 

2. the microphysics of the non-gravitational sector: 
properties of the particles and their interactions 
(mostly the weak interaction) and nuclear physics 
(properties of nuclei and reaction cross-sections), 

3. the matter content of the Universe, mostly through 
the number of baryons per photon 77 = n\j/n^. 

A key physical quantity is the freeze-out temperature of 
the weak interaction, Tf , which can roughly be estimated 
from the Friedmann equation as 

so that we expect the predictions to depend on the num- 
ber of relativistic degrees of freedom, , the Fermi con- 
stant, Gf, the Newton constant G but also the nuclear 
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quantities such as the lifetime of the neutron, the neutron 
to proton mass difference, as well as the fine structure 
constant. The subsequent nucleosynthesis is delayed be- 
cause of the small deuterium binding energy and depends 
on the cross-sections of a dozen nuclear reactions. 

As a consequence, BBN has been extensively used to 
constrain various derivations from the standard physical 
models such as scalar-tensor theories [SHI], variation of 
fundamental constants [SHE] , the existence of extra rela- 
tivistic degrees of freedom, or the existence of particles 
decaying during BBN (see e.g. Ref. [S] for a review on 
BBN and physics beyond the standard model). We refer 
to Refs. [ini E] for reviews on the basics of BBN and on 
the constraints it sets on cosmology and physics. 

Many of these extensions have been motivated by a 
long lasting problem of compatibility between ^Li theo- 
retically predicted abundances and observations. When 
using for 77 the value determined by the cosmic microwave 
observations, the BBN predictions for ^He, D and '^He 
are in very good agreement with those deduced from ob- 
servations. However, there remains, a yet unexplained, 
discrepancy of a factor wS, between the calculated and 
observed ^Li abundances, that has not been reduced, nei- 
ther by recent nuclear physics experiments, nor by new 
observations (see § In] below for a detailed discussion of 
the lithium probleni)T^ 

There are two main channels to produce ^Li. The 
first channel, ^He -I- ^ ^Li + 7, is dominant for 
?7 < 3.3 X 10^1" while the second, '^He + ^Hc '^Bc + 7 
followed by a /3 decay with a half-life of 53 days, is 
dominant for 77 > 3.3 x 10~^° (hence for WMAP value). 
A possibility to modify the abundance of ^Li for the 
larger values of rj is to inject neutrons during the late 
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stages of BBN [SI [T^l US] • That would reduce the amount 
of produced ""Be, and thus of ^Li, since it increases its 
destruction due to a more efficient neutron capture. 

Injecting neutrons is indeed not something easily per- 
formed without modifying the laws of nature and in par- 
ticular the matter sector by including physics beyond the 
standard model of particle physics. Interestingly, such an 
idea can be realized by introducing a mirror sector. Such 
a sector is constructed by assuming that the gauge group 
G of the matter sector is doubled to the product G x G' . 
Imposing a mirror parity under the exchange G G' im- 
plies that the Lagrangian of the two sectors, ordinary and 
mirror, are identical so that they have the same particles 
content such that ordinary (resp. mirroiQ) matter fields 
belonging to G (resp. G') are singlets of G' (resp. G). 
They also have the same fundamental constants (gauge 
and Yukawa couplings, Higgs vev). The latter point im- 
plies that the microphysics (and in particular the nuclear 
sector) is identical in both sectors. The two sectors are 
coupled through gravity, and can eventually interact via 
some couplings so that the general form of the matter 
Lagrangian is 

C = Caie, u, d, (/),...) + Ccie' , u\ d' , (f>\ . . .) + £mix- 

Such a sector was initially proposed by Li and Yang [T3] 
in an attempt to restore global parity symmetry and 
was then widely investigated |151 116j . Any neutral or- 
dinary particle, fundamental or composite, can be cou- 
pled to its mirror partner hence leading to the possibility 
of oscillation between ordinary and mirror particles |18j. 
For instance a mixing term of the form £inix c>c F'^^F^'^ 
will induce a photon-mirror photon oscillation, ordinary 
neutrinos can mixed with mirror neutrinos and oscillate 
in sterile neutrinos |16j . Among all the possible mix- 
ing terms, special attention has been drawn [15] to the 
mixing induced between neutrons and mirror neutrons. 
Such a possibility is open as soon as >Cmix contains a term 
oc {udd){u'd'd') + {qqd){q'q'd')] see e.g. Ref.[l9] for de- 
tails. It was also pointed out [HI that a neutron-mirror 
neutron oscillation could be considerably faster than neu- 
tron decay, which would have interesting experimental 
and astrophysical implications. 

This has motivated experimental searches for nn'- 
oscillations [23, which provide the constraint [24l [26] 

Tosc > 414 s 

at 90% C.L. on the oscillation time. Recently, it has been 
improved :25] to Tosc > 448s ( 90% C.L). 

From a cosmological point of view, mirror particles 
have been advocated as a dark matter candidate 



In particular, mirror baryons do not interact with pho- 
tons and have the same mass as ordinary baryons. From 
our world they can thus be considered as stable, self- 
interacting dark matter particles. The interest in this 
candidate has been revived after it has been argued that 
mirror photons could provide an explanation of the direct 
detection experiments |22| . 

This demonstrates that such a sector is well-motivated 
from theoretical, experimental and cosmological points 
of view. 

Coming back to our primary interest, namely BBN, 
the mirror electrons, photons and neutrinos act as 
extra-degrees of freedom. It is thus clear that it will 
modify the abundance of the light elements in the 
ordinary world and also that one needs to compute the 
abundances in the mirror world. In this framework, 
the abundances of the light elements are governed by 
4 parameters {?7, 77', Tosc, a^}, namely the two photon to 
baryon ratios, 77 and 77', the ratio between the photon 
temperature of the two worlds today, x, and the oscilla- 
tion time Tosc- Two main effects have to be taken into 
account: (1) the fact that the mirror sector accounts for 
extra-relativistic degrees of freedom, hence modifying 
the expansion rate during BBN, and (2) the oscillation 
between neutrons and mirror neutrons. Note again that 
the nuclear physics in the two worlds is not modified 
and strictly identical. 

Early investigations have mostly relied on the first 
effect and studied the nucleosynthesis in the mirror 
world |21j . If the temperatures in both worlds are identi- 
cal, then the effect of the mirror world is equivalent to an 
effective number of neutrino families, SN^s = 6.14 [27) . 
too large a number to be compatible with observation. 
The goal of this article is thus two-fold. First, we want 
to investigate the effect of a mirror sector on the BBN 
predictions. In particular, we will show that it allows 
to set constraints on the parameter Tqsc that are com- 
plementary to laboratory experiments. Then, secondly, 
within these limits, we shall show that there exists a 
range of parameters that alleviate, and eventually solve, 
the ^Li problem. 

The article is organized as follows. We start by recall- 
ing the observational landscape that sustains the lithium 
problem in § [IT] and then describe the implementation of 
a mirror sector in BBN computation in § III The effect 
of the mirror neutron migration is investigated in § IV 
and the predictions are then compared with observations 
in§|Vl 



II. THE LITHIUM PROBLEM 



Any field {j>, fundamental or composite, in the ordinary sector is 
associated to a field <j)' in the mirror sector so that n' corresponds 
to mirror neutrons, p' to mirror protons etc. 



Historically, ^Li obtained a cosmological status with 
the discovery of the Spite plateau [29]. From an obser- 
vational point of view, its abundance keeps a value con- 
stant over a large range of metallicity, (and effective sur- 
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face temperature,) see e.g. Refs. [5DH57] . hence taken as 
the primordial value. From a theoretical point of view, 
as long as one sticks to standard BBN, the computa- 
tion of the abundance of ''Li requires only the knowledge 
of the baryon density, that is nowadays determined to 
an unprecedented precision from analyses of the CMB 
anisotropics and is given, based on the WMAP 7-year 
analysis PSI (see also the WMAP 9-year analysis |55] ) 
by, 

r!b-/i' = (2.249i°;°57) x 10"' (2-1) 
corresponding to: 

V = WMAP = {G-lGtolD X 10-1°. (2.2) 

With this value, standard BBN gives an abundance 
^Li/H of (5.24 ± 0.5) x 10-^° ^ [77|, which is 
considerably higher than almost all observational de- 
terminations. The value determined in Ref. [40 is 
Li/H = {1.23toit) X 10-1° and similarly, the recent 
analysis of Ref. gave Li/H = (1.58 ± 0.31) x 10-^°. 
This leads to a clear discrepancy between theory and 
observation, the stellar observations leading to a value 
too low by a factor of 3-5. 

Many studies have been devoted to the resolution of 
the so-called ^Li problem and many possible "solutions" , 
none of which fully satisfactory, have been proposed. 

First, the hypothesis underlying the standard compu- 
tation and the stellar processes have all been reexamined 
carefully: {i) the nuclear reaction rates used in BBN cal- 
culations have been reexamined [lTtt45j : (ii) the possi- 
ble depletion of ''Li in stars has been discussed [35Hi5] : 
{iii) the temperature scale used in the ''Li abundance 
determination has also been the subject of a lively de- 
bate ISnjHH]; and {iv) the 3-dimensional stellar modeling 
and non-local thermodynamical equilibrium effects have 
also been investigated but the results are not very dif- 
ferent from those determined using 1-dimension models 
with local thermodynamical equilibrium 50 . 

Another possibility is that the solution of the ^Li 
problem could require new physics beyond the standard 
model. Many models have been developed, among which 
{i) the decay of a massive particle during or after BBN 
that could affect the light element abundances and 
potentially lower the ''Li abundance [TH [T31 . Re- 
cently, (ii) an axion condensate which cools the photon 
background leading to value of rj smaller than t^wmap at 
the time of BBN has also been advocated [62l [63] . Exotic 
solutions involving the possibility of a variation of the 
fundamental constants has also been debated 171 IM]. 
Note also the possibility of a modulation of the BBN 
conditions between, (?) the local Universe, where the 
stellar abundances are measured, and (ii) the large scale, 
on which the CMB is measured, so that the lithium 
problem can reflect a breakdown of the Copernican 
principle [65l E^- For a complete analysis, see the 
proceedings of the meeting "Lithium in the cosmos" |67j . 




FIG. 1. Abundances of ''He D, ^He and '^Li as a function of 
the baryon over photon ratio rj (blue) [77] . The vertical stripe 
corresponds to the WMAP baryonic density while the hori- 
zontal area represent the adopted observational abundances 
(green region); see text. The dot-dashed lines correspond to 
the extreme values of the effective neutrino families compati- 
ble with *He observations, N^a = (2.89,4.22). 



Despite this problem, standard BBN has been ex- 
tremely successful and both ''He and D/H show a good 
agreement between the theoretical predictions and their 
observational determinations (Fig. [T]) . 

At T] = T/wMAP, the *He mass fraction is 1^ = 0.2476± 
0.0004 [39l[77]. On the observational side, the determina- 
tion of the helium abundance in extragalactic Hii regions 
is somewhat difficult [S5], Aver et al. [70] find a value 
with large uncertainties, Yp = 0.2534 ± 0.0083, consis- 
tent with the BBN prediction. Using these observational 
limits, one can even set limits on the effective number of 
neutrino families (as defined in the next section) to be 
2.89 < iVeff < 4.22 (red, dash-dotted curves in Fig.[l]). 

The observed deuterium abundance is also in a rea- 
sonable agreement with the BBN prediction which, for 
ri = 77WMAP, is D/H = (2.59 ± 0.15) x 10-^ [SS |77]. 
The weighted mean abundance of deuterium is D/H = 
(3.02 ± 0.23) X 10-^ deduced from the observation of 
about 10 quasar absorption systems (for more details 
see Ref. [71]). Note that the individual measurements 
of D/H show considerable scatter and it is likely that 
systematic errors dominate the uncertainties. Moreover, 
the predicted abundance is somewhat lower than the ob- 
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served mean one. Most of the measurements available 
in the literature have been reviewed in Ref. |72j and 
more recently in Ref. [75]. However, the observation of 
a damped Lyman-a at Zabs = 3.049 has permitted [74] 
a new determination of the abundance of deuterium 
D/H = (2.535 ±0.05) x 10"^ leading to a mean determi- 
nation lower than the previous one [(2.60±0.12) x 10~^]. 
However, since the Hi Ly-a absorption associated to this 
system is redshifted exactly on top of the Ly-a-Nv blend 
emission from the quasar, the errors on this measurement 
arc probably underestimated. Indeed the emission by the 
quasar cannot be recovered easily and is definitively de- 
generated with the exact profile of the absorption. Tak- 
ing into account this observation, Cyburt et al. [5^ find 
a mean value of (3.01 ±0.21) x 10~^ which is comparable 
to the one taken in our calculations. 

However, it is important to note that the high disper- 
sion of the observed D/H values at z '-^ 3 (corresponding 
to t = 3—4 Gyr) spans typically the (2.4— 4.) x 10"^ range. 
Different star formation histories in the galaxies associ- 
ated with the DLAs can explain this dispersion. D is a 
very fragile isotope (it is destroyed in stars at T ^ 10^ K) 
and its destruction rate is dependent on the gas to total 
mass of a galaxy ratio as shown in Ref. [HI] (see their 
Figure 3). As a consequence, D/H destruction factors 
can range between 2 to 10. So, it is reasonable to con- 
sider that these pristine structures are in different evo- 
lution stages and consequently have different gas masses 
and, in situ, processing histories. Finally, as D can only 
be destroyed in the course of chemical evolution, it is 
reasonable to relate the highest D/H observation to the 
primordial post-BBN abundance. Accordingly, we con- 
sider that the primordial D/H abundance can be higher 
than the averaged value and reasonably reach a value of 
4.0 X 10~^ (green, dash-dotted line in Fig. [l|. For a re- 
cent analysis of the deuterium observations, we refer to 
Ref. [71] ■ 



III. RADIATION ERA AND BBN WITH 
MIRROR MATTER 

A. Evolution of the temperature of the mirror 
world 

The Friedmann equation contains both the ordinary 
and mirror matter so that during the radiation era it 
takes the form 

H'=^-^ip. + p'.) (3.1) 

where G = is the Newton constant, p,- the radia- 

tion energy density and where the contribution p'^. of the 
mirror sector has been included. H = d/a is the Hubble 
expansion rate of the scale factor a, a dot referring to a 
derivative with respect to the cosmic time. The energy 
density of the radiation can be expressed in terms of the 



temperature T of the photon bath as 

Pr(T) = ^5*(r)r^ (3.2) 

where (T) is the effective number of relativistic degrees 
of freedom at the temperature T [T]. The determination 
of g*{T) follows the standard computation and has to 
include all relativistic particles. We recall that for any 
species of mass rrii , its energy density is given by 

with + for fermions and — for bosons, and with zi — 
rui/T and u = E/T. As depicted on Fig. [2j two ef- 
fects are important at energy scales relevant for BBN: 
(1) the decoupling of the neutrinos after the freeze-out 
of the weak interaction, which implies that only photons 
are reheated during the electron-positron annihilation so 
that the temperature of the neutrinos differs from T. 
At temperatures smaller than the electron mass me they 
are related by the relation 

after the freeze-out , and (2) the particle-antiparticle 
annihilation (including electron-positron and muon- 
antimuon). 








Log[T/l MeV] 

FIG. 2. The evolution of and g» as a function of temper- 
ature T. At high temperature these two quantities are equal 
but differ after the neutrinos have decoupled. The second 
bumpy drop is due to the muon-antimuon annihilation at a 
temperature of order me. 

Identical relations indeed hold for the mirror world. It 
is clear that if T^^ = Tq today then T' = T at all time 
so that the effect of the mirror world would be to double 
the numbers of relativistic of freedom, which would be in 
obvious disagreement with BBN. It follows that in order 
to be compatible with BBN observations, one needs T' < 
T (i.e. p[ < p^) We thus set 

= xTo (3.5) 
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where a; is a new parameter of the model, which cannot be 
expressed in terms of previous ones. As long as T, T' <C 
trie the temperatures in the two worlds have the same 
scaling with redshift so that T' — xT. However because 
the electron-positron annihilation and the decoupling of 
the weak interaction do not happen simultaneously in the 
two worlds, we expect T' ^ xT in general. We thus set 



T'[x,T]^a^{T)xT, 



(3.6) 



with ax{T) = 1 for T <C m,,. 

Note that because the neutrinos are not reheated dur- 
ing the electron-positron annihilation, their temperature 
always scales as (1 + z) in both worlds so that 



T!. = xT„ 



(3.7) 



at all times. The relation between T' and T is thus ob- 
tained by solving 



Tl{T') = xT,{T). 



(3.8) 



Figure [3] compares the evolution ofT^/T' as a function of 
T to T^T for X < 1. Since the decoupling occurs before 
in the mirror world, the curves are shifted to higher T. 
If plotted as the function of instead of T, the curves 
will be identical up to a translation. 




Log[T/l MeV] 

FIG. 3. T^/T' as a function of T to T^/T for x = 1. (black, 
solid) X = 0.5 (red, dashed) and a; = 0.1 (purple, dotted). 

This translates to the evolution of ax{T) as shown 
on Fig. |4j We recover the limiting behaviours that we 
have discussed above. The deviation from ax — ^ occurs 
when one of the two temperatures is of order of me and 
reaches a typical magnitude of (20-30)%. The amplitude 
increases for smaller x but then the contribution to the 
energy density is suppressed by a factor a;^ so that the 
effect on the Friedmann equation becomes negligible. 

To finish, note however, and contrary to what has been 
assumed e.g. in Ref. [21 that even though at tempera- 
tures relevant for BBN, the ordinary and mirror neutrons 
are non relativistic and do not interact with photons and 
neutrinos, the entropy, s = '^q*T^, of the two worlds 
are separately conserved but that s'/s does not remain 
constant. Indeed, 



s 
s 



al{T)x\ 



(3.9) 



Log[T/l MeV] 



FIG. 4. The evolution of the coefficient ax{T), defined in 



Eq. (|3^ for X = 1 (black), x = 0.8,0.4,0.1 (red dashed, dot- 
dashed, dotted) and x = 1.2,2 (purple, dashed, dotted). As 
discussed above, (T) = 1 at low and high temperatures. 



B. Cosmic dynamics in the ordinary world 



The Friedmann equation (3.11 takes the form 

2 



9 Stt^ , , r 



K{T,T')x^] 



where 



Pr(T) g.^{T) 



(3.10) 



(3.11) 



is a factor that takes into account that a priori T' 7^ 
T. The mirror sector appears as extra relativistic de- 
grees of freedom. Seen from the ordinary world, this 
amount of extra-radiation is characterized by Kx{T) = 
K[T,T'{x,T)]. The deviation from the standard Fried- 
mann equation in the radiation era is often parameterized 
in terms of an effective number of neutrino species, N^e 
as 



Pr 



1 
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off 
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7r2 r4 



(3.12) 



from which SN^tf is constrained to 2.89 < N^s < 4.22 
(§|ll|. As depicted on Fig. [5] (left), this implies that the 
effective number of neutrinos is temperature dependent 
and given by 



5N, 



cff 



4 g,{T) 



Kx{T)x^ 



(3.13) 



7 {TJTY 

At low temperature, it reduces to SN^ii — 7.39i4ra;^, while 
at high temperature it gives SNcH — Q.lAKx'^, as used 
in Ref. [JT]. For T <C mg, it reduces to the standard 
parameterization 

4/3 

1 + - I — I iVoff 



Pr = 



P7 



from which the analysis of the WMAP data combined 
with baryonic acoustic oscillations sets [38] 2.8 < N^e < 
5.9 at 95% C.L., compatible with the standard prediction 
iVoff ~ 3.046. 
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It is clear that if T' = T then, because the two matter sectors are identical, x = 1 and K{T, T') = 1 so that 
SNcff ~ 7.39 or 6Nes ~ 6.14 at low and high temperatures respectively. This is by far incompatible with the existing 
constraints mentioned above. As long as x is too large, the effective radiation content of the Universe is also too large 
and BBN will not be successful. As seen from Fig. [5] (right), we need typically that x < 0.6 in order to have a chance 
to be compatible with current observations. 



Z 



-2 -1 

Log[T/l MeV] 




FIG. 5. Evolution of SN^s (left) with temperature for x =1 (black), 0.8,0.6,0.5,0.4 (red), 0.3,0.2,0.1 (purple). The bump 
corresponds to the electron-positron annihilation while the small increase at thigh temperature is due to the muon-antimuon 
annihilation. Evolution of SNcb (right) with x for T <^ me (black) and rrie <^ T <^ (red). The shaded region is 
experimentally excluded, so that we expect x < 0.6. 
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C. Implementation of the mirror sector for BBN 
computation 

As explained in the introduction, the BBN computa- 
tions depend on 4 parameters. The standard parameter 
is the number of baryons per photon, rj. It has to be com- 
plemented by the constant ratio x, defined in Eq. (3.9) 



and the number of baryons per photon in the mirror sec- 
tor, 77', defined by 



, V 3 

Ph = Ph—x 



(3.14) 



where the factor x^ accounts for the difference in photon 
number density. 

Since the physics in the two sectors is identical, the re- 
action rates are the same in both worlds and BBN is the 
same except for the temperature and the baryonic densi- 
ties. In order to allow for neutron oscillations, we include 
in the nuclear network normal and mirror isotopes with 
initial n' and p' abundances given by: 



Y , , — Y !Lrr3 



Accordingly, the mass fractions are such that 



V 3 

— X 



(3.15) 



(3.16) 



is conserved. 

As discussed in the previous section, instead of the ap- 
proximation T' = xT for the temperature in the mirror 



sector and (5iVoff = 6.14x* to account for the acceler- 
ated expansion, we followed mirror electron and positron 
annihilations to calculate T' [Eq. (3.6)] and added the 
corresponding mirror radiation density, p', in the Fried- 
mann equation [Eqs. (3.1 1 and ( |3.10 )]. This does not 
change the qualitative global features of the results, as 
compared to previous analyses, but leads to significant 
quantitative effects. 

The network includes 16 isotopes and 27 reactions. 
These are the usual 8 isotopes and their mirror partners, 
the 13 main BBN reactions and their mirror counterparts 
plus the n n' oscillation term, described in the next 
section. The reaction rates are the same for ordinary and 
mirror reactions and are those used in previous works e.g. 
Ref. [77], i.e. mostly from Ref. [75], except for ^II(n,7)D 
[7S] and 3He(a,7)^Be [7S]. 

As an example, Fig. [6] shows the result of a BBN cal- 
culation with rj = ?7wMAP, = 10~^° and x = 0.2. 
Nucleosynthesis proceeds in both worlds and the abun- 
dances in the mirror world are consistent with the early 
estimations of Ref. Jl] (e.g. their Fig. 3). Since they 
are of no interest in this study (they do not modify the 
ordinary abundances and are not observable), the abun- 
dances of mirror isotopes other than n', p' and ^He' are 
calculated but not displayed. 

The abundances of the mirror isotopes are globally re- 
duced [Eq. (3.15)] with respect to the normal ones. How- 
ever, it can be seen that after a time of order 300 s, 
the abundance of the mirror neutrons, n', remains much 
higher than the normal neutron abundance. This is due 
to the fact that, for this choice of parameters (T' and p^^), 
mirror BBN is limited to n' decay while normal BBN 
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BBN in both sides (rj =WMAP, r\ '=1 and x=0.2) 
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FIG. 6. Standard BBN (solid lines) and mirror BBN (dashed 
lines, only n , p' and ^He' are depicted, respectively in black, 
red and blue). This computation assumes that the two worlds 
do not communicate, i.e. that Tosc ~ oo. 



is in its full development. It is well known [SI [TH [T3] 
that late time neutron injection in BBN would allevi- 
ate the lithium problem by ^Be destruction through the 
^Be(n,p)^Li(p,a)^He chain. At WMAP baryonic density, 
''Li is produced indirectly by ■^He(Q;, 7) ''Be, that will de- 
cay much later to ^Li. This could be achieved, if some 
fraction of the mirror neutrons would migrate out of the 
mirror world into ours, as we shall now investigate. 



IV. NEUTRON OSCILLATIONS 

In standard BBN, the weak interaction maintains the 
thermal equilibrium between neutrons and protons until 
their rates become slower than the Hubble expansion rate 
at a typical temperature of T « 3.3 GK [79]. This freeze- 
out of the weak interaction is followed by neutron free 
decay until T « 0.9 GK when nucleosynthesis begins. 
This is during this phase of free decay that a mirror world 
allows for neutron oscillation between the two worlds. 
This interval of temperatures would be even wider in the 
mirror world if, as considered here, T' < T and < ph 
(see Fig. [6]). 

To introduce neutron oscillations, it is assumed that 



the mass matrix is given in the form 
m ~ i /2t 1 / Tqj 



M = 



I/tosc m-i/2T 



(4.1) 



where r is the neutron lifetime. It follows that the mirror 
neutron abundance evolves as e"*/'' cos^(i/rosc) [EH] so 
that 



An'^n = - — log COS 

dt V T„„, 



2 / t 
tan I — 



(4.2) 



This rate looks like a decay rate Tosc/2 but modu- 
lated by tan(t/Tosc)- It diverges for t/Tosc = 7''/2 
but (putting aside numerical details), keeps the product 
ATn' (t) tan (t/Tosc) finite, i.e. 

(4.3) 



hm Xn' (t) tan {t/ Tosc) = 

t/To=c->7r/2 



if there is no other source of mirror neutrons. This is 
indeed the case during free mirror neutron decay as we 
neglect An->.n'- When Xn{t)^Xn'{t), at small t, An-i-n' 
is suppressed by a factor tan(i/rosc)- Later, when this 
factor is no longer negligible, Xn{t)^Xn'{t) so that we 
can neglect neutron diffusion to the mirror sector (but 
not the reverse). 



BBN in both sides (r|j„=WMAP, Tij„'=l and x=0.2) 
10^ lo' 




FIG. 7. Same as Fig.lGjbut allowing neutron oscillations with 
Tosc = 3 X 10'^ s: the Be-f ^Li abundance is significantly re- 
duced. 

Fig. [7] shows the result of a calculation with the same 
parameters as in Fig. [6] except that we assume oscilla- 
tions from n' to n with Tosc = 3 x 10'^ s in Eq. (4.2), a 
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value well above the experimental lower limit. Initially, 
because of the difference in baryonic densities, the neu- 
tron abundance is much higher in our world but, dur- 
ing standard BBN the neutron abundance decreases very 
rapidly compared to mirror BBN because of the higher 
temperature and density. Then, when the n abundance 
dropped below the n' one, n injection from oscillating 
n', at the time of ^Be formation, leads to its destruc- 
tion due to a more efficient neutron capture. On the 
contrary, the abundance of D increases: we obtain a re- 
duction of the ''Li -I- ''Be abundance at the expense of a 
moderate higher deuterium. We can observe a fast de- 
crease of the mirror neutron abundance, down to « 0, 
around t « Tosc, followed by a small increase. This late 
time effect is due to the nucleosynthetic n' production by 



the 2H'(d',p')^H'(d',n')'*He' chain when the oscillation is 
turned off. This should not affect the results because it 
happens when the ^Li -f- ''Be abundance has reached a 
plateau and because, in any case, these are fresh mirror 
neutrons for which the oscillation clock has to be reset. 

Note that, before freeze-out the weak interactions, 
that equilibrate the neutron and proton abundances, pre- 
vent any oscillation. We have assumed here that after 
freeze-out, and until p'(n',7')d' reaction becomes effi- 
cient, neutrons are free to decay and oscillate. 



V. COMPATIBILITY WITH OBSERVATIONS 



In the previous section, we have shown that the ''Li primordial abundance can be reduced significantly for values 
of the parameters rj' and Tqsc that were found to be representative of the favorable solutions. We now explore this 
parameters space and compare the predictions to the BBN for ^He, D, ■^He and ^Li primordial lithium (see §0- 




0.0 0.2 0.4 0.6 0.8 1.0 2000 4000 6000 8000 10000 0.0 0.2 0.4 0.6 0.8 1.0 



Logh'io] T„„ (s) Log[^'lo] 

FIG. 8. (left) Limits on 77^9 1/x provided by observational constraints on D (solid; the green is acceptable, i.e. 2.79 x 10~^ < 
D/H < 3.25 X 10"^) "He (dashed; the yellow region is acceptable, i.e. Yp = 0.2534 ± 0.0083) and ''Li (dotted; the light blue 
strip is acceptable, i.e. Li/H — (1.58 ± 0.31) x lO"'^"), allowing for neutron oscillations with Tosc = 3 x 10^ s constant. The 
parameters space in the blue strip reconciled both 'Li and "He, but not D. This conclusion is changed if one allows for a larger 
value of D/H. The thin black lines corresponds, from right to left, to D/H = (3.8, 4.0, 4.2, 4.4, 4.6) x 10"^, which shows that for 
D/H e [3.8,4.6] X 10"^ there exists a region of the parameters (a;,?7io) reconciling 'Li and "He and D. (middle) Same but as a 
function of Tosc and 1/x, keeping jyio = 3 constant, (right) Same but as a function of 77^0 and Tosc, keeping 1/x = 3 constant. 



In 

of Tn 



Fig. Is] (left), we first explore the parameters space 
^ < rf < 10~^ and 0.2 < a; < 1 for a fixed value 
: 3 X 10^ s, assuming neutron oscillations as in 
It shows that for a wide range of parameters, ''Li 



4.2 



Eq, 

and^^e abundances are compatible with observations 
(overlap between the blue and yellow areas), but at the 
expense of D (the green area) . [The constant 1 /x lower 
limit given by ''He just reflects the SNcs < 1.22 limit 
(from § |ll|]. However, it also shows that, if we rather 
adopt a primordial D/H abundance of sa4xl0~^, which 
is reasonable as discussed in § [iTj there exists a region in 
the parameter space {x,ri') for which one obtains a full 
compatibility for all, ''He, D and ''Li isotopes. This is 
confirmed by the Figs, [s] (middle and right) leading to 



the conclusion that there exists a region in the param- 
eters space (a;, Tosc) ioi which the lithium- 7 problem 
is solved, at the expense of a moderately higher D/H. 
This is, in fact, a common feature that models succeed- 
ing in reducing the ''Li BBN production come with an 
over production of D. This nevertheless can be cured e.g. 
by non-thermal photons injection by the radiative decay 
of long-lived particles after the BBN g^. Fig. [s] (right), 
shows that if we adopt the mean value observed today, 
D/H = (3.02 ± 0.23) x 10-^ for any Tosc values between 
10^ and 10^ s there are combinations of rj' and x, that 
reconcile D/H and helium-4 (the green area). 

To finish, we let the oscillation time vary and explore 
the excursion of D/H allowed when rj' is varying between 
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10~^° and 10~^ and x between 0.2 and 1, while 77 is fixed 
to its WMAP value. We then determine the range of 
values of D/H that would be required in order to have 
at least 1 model (i.e. a set of x and 77') so that D/H, 
lithium-7 and helium-4 are in agreement with their spec- 
troscopic observations. Figure |9] shows that such a solu- 
tion exists only if Tosc < 6.7 x 10^ s and would require 
a typical value of D/H of the order of 4 x 10~^. For 
larger oscillation times, the mechanism is no more effi- 
cient and plays no role on BEN. It is important to note 
that with the current status on the observation of D/H 
(i.e. D/H = (3.02±0.23) x 10"^), no global reconciliation 
is allowed, whatever the value of the oscillation time. 



l<r|j„'<10 and 0.2<x<l. 



1 2 3 456789 10 

\tz (10^ S) 

FIG. 9. As a function of the oscillation time Tosc, we let x vary 
between 0.2 and 1 and rj' between 10"^" and 10~^, and explore 
the value of D /H that can be reached while rj — t^wmap (blue 
dashed line). Within this range, we find models, i.e. at least 
one set of parameters (77', a;) such that D/H, '^Li and Yp would 
agree with their spectroscopic abundances (black line). This 
region is a narrow strip around D/H = 4 x 10~^ (dot-dashed 
line) and exists only if Tosc < 6 x lO'^ s. The green strip 
corresponds to D/H = (3.02 ± 0.23) x 10~^ a range for which 
the model does not allow for any resolution of the lithium-7 
problem, whatever the oscillation time. 



the light elements synthesized during BBN. The main 
motivation lies in the fact that such a model allows for 
an injection of neutrons during the late stages of BBN, 
an effect that in principle shall lead to a decrease of the 
abundance of lithium-7. It is thus a good candidate in 
order to resolve the lithium problem, i.e. the fact that the 
abundance of lithium-7 predicted under the hypothesis 
that the baryon density is given by the CMB observations 
is about 3 times larger that the value obtained from the 
spectroscopic observations. 

We have implemented mirror neutrons and neutron os- 
cillation in a BBN code. The model depends on 3 extra- 
parameters, X the relative mirror photons temperature 
relative to the CMB temperature today, Tqsc the oscil- 
lation time and 77' the mirror baryon number. We have 
shown that as soon as a; < 0.6 the helium-4 abundance is 
in agreement with the observations while the agreement 
for both the deuterium and helium-4 abundances can be 
obtained within a large domain of the parameters space. 
This is a new example of the power of BBN to constrain 
deviations of the standard model of particle physics. 

Concerning the lithium-7 problem, we have shown that 
if D/H = (3.02 ± 0.23) x lO'^ then D/H, helium-4 
and Li/H cannot be reconciled whatever the parameters 
(x, 7?', Tosc)- However, if D/H e [3.8,4.6] x 10"^ then 
the abundances of all light element (helium-4, deuterium 
and lithium-7) can be reconciled for a 77 = 77vvmap- In 
conclusion, there exists a region in the parameter space 
(x, 77', Tosc) for which the lithium-7 problem is resolved, 
at the expense of a moderately higher D/H. Indeed such 
a value is not unreasonable, as discussed above. Note 
that this requires that Tosc < 6.7 x 10^ s. This upper 
bound on the oscillation time is 15 times larger than the 
lower bound provided by laboratory experiments so that 
an other way to falsify this explanation would be to close 
the window on Tqsc- 

In conclusion, mirror neutrons do alleviate the lithium- 
7 problem and neutron oscillation offers an interesting 
class of falsifiable models that can potentially solve the 
lithium problem, an issue that can be resolved once spec- 
troscopic measurement of D/H are available or better 
lower bounds on the oscillation time are established. 
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